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Abstract 



It has been suggested by Groot Nibbelink and Pospelov that Lorentz invariance 
can be an emergent symmetry of low-energy physics provided the theory enjoys a non- 
relativistic version of super symmetry. We construct a model that realizes the latter 
' symmetry dynamically: it breaks Lorentz invariance but leaves the supersymmetry 

' generators intact. The model is a super symmetric extension of the dynamical aether 



theory of Jacobson and Mattingly. It exhibits rich dynamics and possesses a family of 
inequivalent vacua realizing different symmetry breaking patterns. In particular, we 
find stable vacua that break spontaneously spatial isotropy. Supersymmetry breaking 
terms give masses to fermionic and bosonic partners of the aether field. We comment 
on the coupling of the model to supergravity and on the implications for Hof ava gravity. 



1 Introduction 

Recently the idea that Lorentz invariance (LI) may not be a fundamental symmetry of nature 
has gained considerable attention. Indeed, violation of LI is present in one form or another in 
many theories of gravity devised to solve the problems of Einstein's general relativity (GR). 
This is the case of the infrared modifications of gravity, which attempt to address unresolved 
problems of cosmology such as the nature of dark matter and dark energy. Examples are 
the ghost condensate [1], Lorentz violating massive gravity [2], the Galileons [3] and kinetic 
gravity braiding [4] . It has been shown [5] that violation of LI allows a simple and technically 
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natural explanation of the dark energy. Also breaking of LI can be a consequence of quantum 
gravity. In particular, in the approach to quantum gravity recently proposed by P. Horava 
[6] (see also [7, 8, 9, 10]) one manages to construct a potentially renormalizable theory of 
quantum gravity at the price of abandoning LI at very high energies. 

Of course, LI is one of the best tested symmetries experimentally [11, 12, 13] and one 
may wonder how this can be reconciled with the ideas mentioned above. In other words, in 
order to be viable, any Lorentz breaking model must incorporate a mechanism that ensures 
the emergence of LI to a very high accuracy, at least within the Standard Model (SM) of 
particle physics at relatively low energies. 

In [14, 15] it was suggested that such mechanism can be provided by supersymmetry 
(SUSY). Namely, consider the standard SUSY algebra and remove the boosts from it [16, 17] 
1 . The remaining generators still form a closed algebra that we will call non-relativistic SUSY 
(or, where it will not lead to confusion, just SUSY for short). As will be discussed below, this 
is actually the most general minimal superalgebra containing spatial rotations and space-time 
translations. A theory invariant under the non-relativistic SUSY is in general not Lorentz 
invariant: it is possible to explicitly construct supersymmetric terms in the Lagrangian that 
violate LI. However, it turns out that with the field content of the Minimal Supersymmetric 
Extension of the SM (MSSM) all these Lorentz violating operators are of dimensions five or 
higher 2 . In other words, LI emerges as an accidental symmetry at the renormalizable level. 
The eventual breaking of SUSY introduces Lorentz violation also at this level, but the effect 
is within the existing bounds, provided the masses of the MSSM superpartners lie sufficiently 
below the Lorentz violation scale [14, 15]. 

The analysis of [14, 15] is performed in the flat space-time where violation of LI cor- 
responds to the existence of a globally defined preferred frame. However, when gravity is 
taken into account the space-time becomes dynamical and it is clear that the preferred frame 
must become dynamical as well. Thus to implement the mechanism of [14, 15] in gravita- 
tional models with Lorentz violation it is necessary to construct SUSY extensions of theories 
describing the preferred frame dynamics. 

In the present paper we consider a representative theory of this class - the Einstein-aether 
model [18, 19]. In this model the preferred frame is defined by a vector field u m (aether) 

1 Refs. [14, 15, 16, 17] further restrict the SUSY algebra stripping off from it the spatial rotations as well. 
Here, we prefer to keep the invariance under rotations explicit and thus retain the corresponding generators 
in the symmetry algebra. 

2 Dimension five operators can be forbidden by imposing additionally CPT invariance, leaving dimension 
six terms as the lowest Lorentz violating contributions. 
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which is constrained to be time-like and have unit norm. We present a supersymmeric 
extension of the model and analyze its consequences. We restrict to the case of global SUSY 
deferring the coupling of the model to supergravity (SUGRA) for the future. This is a good 
approximation when the energy scale characterizing Lorentz violation is small compared to 
the Planck mass. Interestingly enough, we will find that the fluctuations of the aether field 
around its vacuum expectation value (VEV) also exhibit LI at the level of low-dimension 
operators. This extends the mechanism of [14, 15] to the supersymmetric aether sector. 

The construction we present in this paper allows to realize the non-relativistic SUSY as 
the residual symmetry left over by the aether VEV in an otherwise super- Poincare invariant 
setup. In this language the mechanism of [14, 15] can be formulated as follows: so long as the 
aether VEV preserves SUSY, the Lorentz breaking does not propagate into the MSSM sector 
at the renormalizable level. The eventual breaking of SUSY, that must be incorporated in 
any realistic model, is unrelated to the dynamics of the aether. It is assumed to come from a 
different source characterized by a lower energy scale. It is worth mentioning two intersting 
recent papers [20, 21] where a general recipe was given to construct supersymmetric versions 
of a wide class of theories including Lorentz breaking models such as ghost condensation. 
However, in the models constructed in this way, the Lorentz violating backgrounds inevitably 
break SUSY at the same time. This results in the absence of a hierarchy of scales between 
SUSY and Lorentz breaking and so the LI of MSSM is not protected in these models. Our 
approach differs in spirit. Our main goal is to construct a Lorentz violating theory with 
unbroken SUSY and study what kind of restriction this requirement imposes on the structure 
of the theory. 

The paper is organized as follows. In Sec. 2 we describe the Einstein-aether model. In 
Sec. 3 we introduce the non-relativistic SUSY algebra and briefly review the arguments of 
[14, 15]. The SUSY extension of the aether model is constructed in Sec. 4. In Sec. 5 we 
analyze the effects of SUSY breaking on the super-aether. We conlude in Sec. 6 by discussing 
our results and future directions. 

Throughout the paper we follow the conventions and notations of [22]. In particular, 
the metric signature is (—,+,+,-1-); Lorentz indices are denoted by the letters from the 
middle of the latin alphabet, m, n, ■ ■ ■ — 0, 1, 2, 3; letters from the beginning of the alphabet, 
a,b, . . ., are used for purely spatial indices. We work with two-component spinors, the spinor 
indices are denoted by Greek letters with or without an overdot depending on the chirality. 
The spinor indices are raised and lowered using the two- index antisymmetric tensor e a ^\ 
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ip a = e a ^ip/3 } tp a = EapipP , where e = £21 = 1- We will often omit the spinor indices 
adopting the following convention for the multiplication of spinors: iprj = ip a r] a , ipfj = ipaV a - 



2 Review of Einstein-aether 

The action of the Einstein-aether theory has the form [18, 19] 

S = S GR + S m , (1) 

where 

M 2 f 

s GR = -^ V=gRd A x (2) 
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is the standard Einstein-Hilbert action for gravity and the aether action is given by 
M 2 f / 

S„ = —f I V=g(ciV n u m V n u m + c 2 (V m u m ) 2 + c 3 V n u m V m u n 

-c 4 u r u s V r u m V s u m + \(u m u m + l)^jd 4 x . 



(3) 



Here u m is a dynamical vector field and A is a Lagrange multplier that enforces u m to be 
time-like with unit norm, 

u m u m = -1 . (4) 

The dimensionful parameter sets the scale of violation of LI; the dimensionless constants 
ci, . . . , C4 are the remaining free parameters of the theory 3 . Eq. (3) is the most general action 
for a vector with fixed norm containing up to two derivatives. 

Two comments are in order. The third term in (3) can be cast after integration by parts 
into the same form as the second term plus a non-minimal coupling of the vector to gravity, 

c 3 V n u m V m u n ~ c 3 (V m u m ) 2 + c 3 R mn u m u n , (5) 

where R mn is the Ricci tensor. Next, the energy- momentum tensor of the aether is propor- 
tional to Ma,. Thus if the scale of Lorentz breaking is much lower than the Planck scale, 
Ma, M ph the effect of the aether on gravity is small. In this regime it is self-consistent to 



3 This parameterization is redundant: the change of can be absorbed into redefinition of Cj's. One 
could fix the ambiguity by choosing = M p i as in [19]. Then the scale of Lorentz violation would be 
set by the product M p ia. However, we find convenient to keep the scale of Lorentz breaking as an explicit 
parameter while assuming the constants a to be of order one. 
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consider the aether in external metric, neglecting its back-reaction on the geometry. In this 
paper we concentrate on the aether theory in flat space-time. 

The ground state corresponds to a homogeneous configuration of the aether which by a 
Lorentz transformataion can always be cast into the form 

iCc = (l,0,0,0). (6) 

This configuration breaks LI down to the SO (3) subgroup of spatial rotations. The model 
decsribes three propagating degrees of freedom: one helicity-0 and two helicity-1 modes. 
These have linear dispersion relations 

w = S(j) |k| (7) 

with the velocities 4 [19] 

2 _ Ci + C 2 + C 3 2 _ Ci 

S(0) " c 1 + c 4 ' S « - ^T* ' { ) 

The difference in the propagation speeds of the modes with different helicities clearly mani- 
fests the breakdown of LI. 

Generally we can consider coupling the aether to other sectors of the theory including 
SM. This will lead to a strong violation of LI in these sectors. Let us illustrate this point 
on an example of a scalar field 0, that we take as a toy model for the SM sector. Possible 
interactions of <p with the aether include a dimension-four operator 

^u m u n d m ^d r ^. (9) 

In the background (6) this will modify the dispersion relations of the 0-particles, 

u 2 = s 2 k 2 + m 2 , (10) 



s 2 = — . (11) 



where now the maximum propagation velocity s is not equal to one, 

1 

1 + k 

For a single field this modification can be absorbed by redefinition of the units of time or 
space. However, the real problem arises when we consider several particle species. In general, 
their interactions with the aether will involve different coupling constants leading to different 
maximal propagation velocities. Non-observation of such differences within SM leads to very 
tight constraints on the couplings of the form (9) [13]. As we now discuss, the unwanted 
couplings can be fobidden by SUSY. 



4 To obtain (8) from the expressions of [19] one should take the limit a -C 1 corresponding to the case 
when back-reaction of the aether on gravity can be neglected. 
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3 Lorentz invariance from non-relativistic SUSY 



Consider the non-relativistic symmetry algebra consisting of the generators of 3- dimensional 
translations P a and rotations J a , a = 1,2,3, supplemented by the time translation P . Its 
minimal supersymmetric extension contains two supercharges Q a , a = 1,2, transforming in 
the spinor representation of the 5*0(3) group of rotations, and their Hermitian conjugates 
Q a . The most general commutation relations compatible with the SO (3) symmetry and 
Jacobi identities are 5 , 

{Q a ,Qfs}=2Aa a a(3 P a , (12a) 
{Q a ,Q fS } = -2A*(a a )^Pa, (12b) 
{Qa, Q P } = 2B(a a f a P a - 2C5?P , (12c) 

[Pa, Qa] = [Pa, Q*] = [P , Qa] = [P , Q*] = . (12d) 

Here A is a complex constant, B and C are real, and (<7 a )g are the Pauli matrices 6 . The 
commutators involving the angular momenta J a have the standard form dictated by the 
5*0(3) symmetry and we do not write them explicitly. Note that by choosing the dimension 
of the supercharges to be (length) -1 / 2 the constants A and B are made dimensionless. Then 
C has the dimension of velocity. Existence of this universal velocity encoded in the SUSY 
algebra allows to connect units of space and time and, as we are going to see, eventually leads 
to the emergence of LI at low energies. We will refer to the algebra (12) as "non-relativistic 
SUSY". 

We now show that it is equivalent to the standard 4-dimensional super-Poincare algebra 
with the boosts stripped off. Indeed, it is straightforward to check that by an appropriate 
linear transformation of the supercharges 

Qa ^ Qa = dlQa + ^afiQ^ (13) 

the coefficient A can be set to zero. Then, assuming the generic case when both B and 
C are non-zero, we can set B = C = 1 by rescaling the generators P a , P (i.e., choosing 

5 We consider the case of non-zero C in (12c). If C = there is an extra solution with 

[■Fb, Qa] = CiQ a + CiE a pQ® , 

where C\ is real and C2 — —C\A/B. 

6 We limit our consideration to proper algebras, in the sense that the r.h.s. of the commutation relations 
involves the generators only linearly. As was pointed out in [23, 24], one could consider more general 
constructions allowing the coefficients in (12) to depend on the Casimir operators such as P a P a and P$. We 
do not pursue this route in the present paper. 
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appropriately the units of length and time). Finally, redefining 

Q a H- Q a (14) 

one obtains 

{Q a ,Qp} = 2a™P m , (15a) 

{Qa, Qp} = {Qa, Qp] = [Pm, Qa] = [P m , Qa} = , (15b) 

where a d = — This is nothing but the commutation relations of 4d SUSY. 

ap a ° 

We now briefly review the argument of [14, 15] showing that the symmetry (15) leads 
to LI at the renormalizable level of MSSM. As in the relativistic case, the representations 
of the algebra (15) are classified in terms of the superfields: indeed, the boost generators 
never appear in the superspace construction. The MSSM Lagrangian is constructed out of 
a number of chiral matter $(/) and real gauge Vjy) superfields. It involves integration over 
the Grassmann variables 9 a : f d 2 9 in the superpotential term and J d 2 9d 2 9 in the kinetic 
(Kahler) part. Table 1 summarizes the mass dimensions of the objects that can appear in 
the Lagrangian. 



matter field $(/) 


1 


gauge field V (J ) 





gauge field strength W(j) a 


3/2 


space-time derivative d m 


1 


supercovariant derivatives D a , D a 


1/2 


chiral measure J d 2 9 


1 


super-space measure J d 2 9d 2 9 


2 


Lorentz breaking VEV u m 






Table 1: The mass dimensions of various objects entering the MSSM Lagrangian. 

Let us try to construct from these ingredients a Lorentz violating contribution to the 
Lagrangian. This contribution must be invariant under the remaining symmetries of the 
theory such as gauge invariance. It is convenient to adopt a formally Lorentz covariant 
description using the vector u m with a Lorentz breaking VEV (6) as a compensator 7 . Then a 

7 In Sec 4.2 we will see that the SUSY aether can provide richer patterns of Lorentz breaking. For the 
sake of the argument, we concentrate here on the simplest pattern corresponding to the VEV (6). 
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Lorentz violating term in the Lagrangian will contain a number of Lorentz indices contracted 
with u m . The superfields 3>m, Vu\ are Lorentz scalars. Thus the Lorentz violating term 
must involve at least one derivative or the gauge field strength W<j\ a . By inspection, with 
the account of Table 1, one finds that there are no Lorentz violating terms of dimensions 2 
or 3. The only allowed term of dimension 4, that is not a total derivative, has the form 

Cf v = K{IJ) u m J d 2 9 $ (/) <9 m $ (J) + h.c. , (16) 

where I ^ J. However, in MSSM this term is forbidden by the gauge invariance. Indeed, 
all fields of MSSM are charged under the gauge group and the plain derivative in (16) is not 
gauge invariant 8 . One could try to remedy this by introducing the gauge-covariant derivative 

T? m $ = ~a^D & e- v D a e v ^ . (17) 
However, this expression is no longer chiral, 

bp m ® ± o , 

and cannot be used in the superpotential term. Thus one concludes that in MSSM it is 
impossible to write any Lorentz violating operator of dimension 4 or less. 

To avoid confusion let us stress that the above argument does not imply complete restora- 
tion of LI. One can easily construct higher dimension operators that do violate LI. An ex- 
ample is the dimension 6 contribution 

771 Yl P 

Cfv = Jftj dW V m ®V n <S> (18) 

in the kinetic part. Here M LV is the scale of Lorentz violation. Due to the eventual SUSY 
breaking the Lorentz violation from higher-order operators feeds into the lower- dimensional 
ones [14, 15]. In particular, the operator (18) gives rise to the dimension 4 terms of the form 
(9). However, the coefficients in front of these terms are suppressed by the ratio m 2 SUSY /Ml V) 
where m S usY is the scale of SUSY breaking in MSSM (the scale of the superpartner masses). 
This ratio can be small enough to satisfy the experimental bounds on Lorentz violation. 

8 Note that terms of the form (16) do appear if we include into consideration sterile neutrinos that are 
gauge singlets. 
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4 Super-aether 



4.1 The Lagrangian 

We now turn to the construction of the supersymmetric aether theory. We impose the 
requirement that this theory must be compatible with the mechanism for emergence of LI 
reviewed in the previous section. This determines the choice of the SUSY multiplet to 
embed the aether vector u m . Indeed, for the argument of the previous section to hold it is 
necessary that a constant VEV of u m preserves SUSY. This implies that u m must be the 
lowest component of the multiplet. The simplest multiplet with these properties corresponds 
to a chiral superfield U m , 

D & U rn = . (19) 
This superfield is a Lorentz vector. In components we write 

U m = u»(y) + V28 a r]™(y) + 6 2 G m (y) , (20) 

where 8 2 = 6 a 6 a and y m = x m + i9a m 9. Note that the aether vector u m is now allowed to 
be complex, 

u m = u^ + iu^, (21) 

with -u^, u™ real. The vector-spinor field rj™ in (20) is the aether superpartner and G m is 
the auxiliary component. By analogy with Eq. (4) we impose the constraint 9 

U m U m = -1 . (22) 

Note that U m has vanishing mass dimension. 

Following the ordinary aether theory as the guideline we look for the super aether La- 
grangian in the form 

C = Mlt . (23) 

In other words, we factor out the scale of Lorentz breaking from the Lagrangian. As long as 
we are interested in the low-energy physics we can restrict to the operators of dimensions up 
to 2 in C. We also assume that the aether field is the only source of Lorentz violation, and thus 
the Lagrangian must be a Lorentz scalar. With these restrictions the only possible form of 
the kinetic term is a function of the dimensionless combination U m U m because the superspace 
integration already contributes dimension 2, see Table 1. Turning to the superpotential part, 

9 The alternative choice of the constraint U m U m = —1 does not lead to a consistent theory. 
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all terms without derivatives are trivial due to the constraint (22). Moreover, terms with 
one space-time derivative vanish due to the identity 



U m d n U m = 



that is a consequence of (22). Finally, the terms with more derivatives are of higher di- 
mensions. One concludes that at the level of operators of dimensions less or equal 2 the 
superpotential vanishes. In this way we arrive at the most general super-aether Lagrangian 



= Ml 



j d 2 9d 2 9 f(U m U m ) + (^J d 2 9 A(U m U m + 1) + h.c. 



(24) 



where / is an arbitrary function and we have implemented the constraint (22) by a super- 
potential term with a Lagrange multiplier chiral superfield A. 

One makes an important observation. Besides the usual LI that acts on the coordinates 
and the indices of the field U m , the Lagrangian (24) possesses an extra SO(3, 1) invariance 
that acts only on the U m indices. In fact, if one treats the index 'm' as internal, one 
immediately recognizes in (24) the Lagrangian of the standard supersymmetric non-linear 
sigma-model with internal Lorentz group. Thus one concludes that at the level of the lower- 
dimension operators the super-aether is equivalent to a Lorentzian sigma-model with Kahler 
potential f(U m U m ). The symmetries of (24) contain the product 

Lorentz invariance x internal 5*0(3, 1) . 

Even when the vector u m acquires non-zero VEV the diagonal subgroup of this product 
remains unbroken. This means that the theory for perturbations around the ground state 
effectively continues to be Lorentz invariant. In particular, as we will see explicitly below, all 
perturbations have the same dispersion relation. Thus LI emerges as an accidental symmetry 
due to SUSY even in the Lorentz violating sector! Let us emphasize that this is happening 
despite that U m is neutral with respect to all (internal) gauge symmetries. Rather, the 
emergence of LI in this sector appears as a consequence of the unit-norm condition. 

Let us stress again that this does not mean that LI is completely restored. Once we go 
beyond the low-energy theory it is easy to write higher-dimensional contributions into the 
super-aether Lagrangian where U m is contracted with space-time derivatives, such as the 
operators 

J d 2 9 {d m U m f + h.c , J d 2 9d 2 9 U m U n dJJi d n U l 

multiplied by the appropriate dimensionful coupling constants. This removes the accidental 
internal S0(3, 1) symmetry and thus leads to a genuine breaking of LI by the VEV of u m . 
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Alternatively we can consider coupling of the super-aether to another sector of the theory. If 
this sector contains (at least two) chiral gauge singlets, it is possible to construct a dimension 
4 coupling 

4v = J d 2 e U m <S> {I) d m <S> {J) + h.c. (25) 

that reduces to the Lorentz violating operator (16) in the constant aether background. If 
gauge singlets are absent, as in the case of MSSM, there are always higher- dimensional 
couplings such as, e.g., 

Cfl = Mil J d2dd2 Q U m U n V rn §V n <& . (26) 

This is the generalization of (18) to the case of dynamical aether. Finally, the internal 
SO(3, 1) will be broken by gravity that mixes the space-time and vector indices through 
covariant derivatives. 

It is instructive to write the Lagrangian (24) in components. Using (20) and a similar 
decomposition 

A = l(y) + V29ay) + 2 H(y), (27) 
after integration over the anticommuting variables we obtain, 

C = dbos "I" £ferm j (28) 

where the bosonic and fermionic parts are 

C bos = Ml \-f mn d r u m O r u n + f mn G m G n + [H{u m u m + 1) + 2lu m G m + h.c] ] , (29) 

m 2 r 

" if nr d m fj n v m Vr 



r 

*-ferm 



+f" (- iu m d r u n ffWrr - iu r d m u r fj n a m r] n - u m G n f\ m f] n + ^(fjmfjn)(v m V n ) 



+f" ( - iu m u n u r d s u T r] m a s r] n - u m u n u s G s r] m r] n + u m u n {rT^rfrf)) 

+ \f""u m U n U s U r (fj m fj n )(V S V r ) - IVmV" 1 - 2(tVm)u m + h.C 



(30) 



Here the "effective metric" is 

and primes denote derivatives of the function / with respect to its argument 

w = u m u m . (32) 
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Let us set the fermionic components to zero and concentrate on the bosonic part (29). 
The equations of motion for the fields G m , I imply that these fields vanish, i.e. they are 
non-dynamical. The field H is a Lagrange multiplier that enforces the unit-norm constraint 
on the aether u m . For purely real aether Eq. (29) reduces to the standard aether Lagrangian 
(3) with 

ci = 2/'(-l) , C2 + c 3 = C4 = . (33) 

Due to the relation (5) considerations in flat space-time do not allow to determine the 
coefficients c 2 , c 3 separately but only their sum. 

4.2 Vacua and fluctuations 

We now classify the ground states of the model in the general case of complex aether. The 
unit-norm constraint implies two equations for the real and imaginary parts of u m , 

u^u Rm - u™u lm = -1 , (34a) 
Mr«i™ = . (34b) 

Consider first the case when the vector u R is time-like. By a suitable Lorentz rotation it can 
be aligned with the time axis. Then due to (34b) u[ n has only spatial components and by a 
3d rotation we direct it along the third axis. Taking into account (34a) we obtain a family 
of ground states, 

M mc = (cos a, 0, 0, isina) , (35) 

parameterized by a G [0, 7r/2]. Note that for a ^ the ground state, besides breaking boosts, 
breaks also the spatial isotropy. Indeed, for a ^ 0, 7r/2 the group of invariance of (35) reduces 
to the SO (2) rotations in the (l,2)-plane. For a = tc/2 the unbroken group is enhanced up to 
SO(2, 1). Only for a = one recovers the 3-dimensional SO (3) rotations. Thus we conclude 
that SUSY has introduced a qualitatively new feature into the aether model: existence of 
a family of inequivalent vacua, generally breaking spatial isotropy together with LI. Let us 
stress, however, that breaking of rotational invariance has a very different physical meaning 
from the breaking of boosts. In a sense, the boost invariance is broken explicitly. Indeed, 
there is no state in the theory where it would be recovered 10 . On the other hand, breaking 

10 One may think that the boost invariance can be at least partially restored, as it appears to happen in 
the vacuum (35) at a = ir/2 and in the vacua (36). However, we are going to see that these vacua contain a 
ghost in the spectrum of perturbations and thus are dynamically inaccessible. 
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of the rotation group is truly spontaneous: the symmetry breaking configurations smoothly 
connect to the vacuum with a = where the full 5*0(3) group is restored. 

Apart from (35) the constraints (34) possess a branch of solutions when both and u™ 
are space-like. Due to the orthogonality condition (34b) we can direct these vectors along 
the second and third axes. Then Eq. (34a) yields 

< ac = (0,0,sh/3,ich/3), (36) 

where f3 > 0. These vacua leave unbroken the subgroup SO(l,l) of the Lorentz group 
corresponding to the boosts in the (0,l)-plane. 

To get more insight into the dynamics let us study small perturbations around the vacua. 
We start with the background (35). Writing 

u m = uZ c + v m , (37) 

one finds that the unit-norm constraint relates the components i>o and V3. It is convenient 
to express them in terms of a single field v\\, 

v = —iv\\ sin a , v 3 — v\\ cos a . 

Substituting this into (29) we obtain the quadratic Lagrangian, 



Cbos = -M' K 

where 



(-«,/' + (!- w 2 )f") d m v\\ d m v {l + f'J2 d ™ v « dmf} °\ ' ( 38 ) 



0=1,2 



w = — cos let . (39) 



Note that this Lagrangian is effectively Lorentz invariant if we treat the perturbation com- 
ponents Vi, t>2, v || as scalars. In particular, all modes have the same linear dispersion relation 
with unit velocity. This is the manifestation of the emergent LI discussed above. On the 
other hand, the symmetry between v» and v±, V2 is clearly broken in (38). The Lagrangian 
is free of ghosts or any other pathologies as long as 

/' >0, -wf +(l- w *)f" >0. 

Therefore, the model allows breaking of spatial isotropy (in the vacuum with a 7^ 0) and 
still possesses completely stable spectrum of perturbations. Note that this situation cannot 
be realized in the real-aether model of Sec. 2. In that case considering space-like aether 
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corresponds to switching the sign in the constraint (4). Then it is straightforward to check 
that the spectrum of perturbations contains a ghost. Instead, for the super-aether model 
changing the sign in the constraint (22) does not affect the physics: it simply corresponds to 
multiplying U m by %. We have seen that the model is stable so long as one of the components 
of the vector - its real or imaginary part - is time-like. 

In the purely spatial background (36) the situation is different and, actually, is analogous 
to the case of spatial real aether. Here the unit-norm constraint relates i>2 and t>3, and the 
quadratic Lagrangian reads 

Uos = -K[- f d mv, d m v + f d mVl d m v x + (wf + (w 2 - l)f ) d m v\\ d m v {l ] , (40) 

where now w = ch2/3 and v\\ is introduced as 

v 2 = —iv\\ ch/3 , v 3 = i>|| sh/3 . 

Clearly, one of the modes is always a ghost. Thus the purely space-like aether background 
is pathological and we do not consider it in what follows. 

Finally we analyze the fermionic content of the model. The spinor field £ in the last line 
of (30) is a Lagrange multiplier that enforces the constraint 

u m Vm, = . (41) 

Solving this constraint in the background (35), 

77o 7 = — ir)\\ y sin a , = rj^ cos a , 

we obtain 



•'ferm 



£ fl , ,„ = ( - wf + (1 - <r J )/")/WT'''// +f'Yl d mVaV m Va + h.c] + . . . , (42) 



0=1,2 

where w is given by (39) and dots stand for the terms with higher powers of the fermionic 
fields. This Lagrangian describes three fermions in the Weyl representation with respect to 
the emergent LI that match the three complex scalars of (38). 



5 Breaking of SUSY 

In this section we discuss the effects of SUSY breaking on the super-aether model. The 
breaking of SUSY is conveniently described using spurion superfields whose higher compo- 
nents have non-vanishing VEVs. We will assume the spurions and their non-zero components 
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to be Lorentz scalars, so that they do not source additional violation of LI. Then the SUSY 
breaking Lagrangian of the lowest dimension is, 

C SB = -Ml j d 2 9d 2 9 [S m g (1) (U m U m ) + S m g m (U m U m )] , (43) 

where g(i), g( 2 ) are arbitrary functions and the spurions have the form 

S {1) =m\ 1) 9 2 9 2 , (44a) 
S {2 ) = m {2) (9 2 + 9 2 ) . (44b) 

The parameters m^, m,( 2 ) have the dimension of mass; they set the scale of SUSY breaking 
in the aether sector. The Lagrangian (43) reads in components 



Csb = —M. 



m fi)9(i)(w) + m(2) (g{ 2) {w)G m u m - ]- g" 2) (w)u m u n (r] m r] n ) + h.c. 



(45) 



where w is defined in (32). This must be combined with the supersymmetric part (29), (30). 

It is straightforward to work out the effect of the first term in (45). Consider the 
aether vacuum (35). One observes that the modulus a acquires a potential proportional 
to 9(i) (— cos 2a). If gLA—1) > the modulus is stabilized at a — 0. Thus we conclude 
that SUSY breaking lifts the degeneracy of the aether vacua and, under broad assumptions, 
singles out the purely time-like aether configuration. Note that this configuration is real and 
preserves the spatial isotropy. The imaginary part of the aether perturbations acquires a 
mass of order 777(1) in this vacuum. 

The effect of the contribution proportional to m( 2 ) in (45) is subtler. At first sight, the 
second term in the round brackets seems to be a candidate for the fermionic mass term. 
However, due to the constraint (41), this term actually vanishes to quadratic order in the 
purely real aether background. Nevertheless, the fermions do get masses, though in an 
indirect way. We observe that the m( 2 )-contribution modifies the equations for the non- 
dynamical fields G m , I, so that the field I no longer vanishes. Instead we obtain in the 
vacuum 



U = ^ ■ (46) 

Substituted into the last line of the fermionic Lagrangian (30) this endows i] m with the mass 
of order 777.(2). 

To sum up, we have found that the general result of SUSY breaking is to generate 
masses for the fermions and the imaginary part of the aether. Below the SUSY breaking 
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scale one is left with the ordinary real aether described in Sec. 2. However, the memory of 
the SUSY origin of the theory is not totally lost: it is encoded in the relations (33) between 
the coefficients of the aether action 11 . 

6 Discussion and outlook 

In this paper we have considered a prototypical model of Lorentz violation with a dynamical 
preferred frame - the Einstein-aether theory - and have constructed the supersymmetric 
extension of the aether sector. Our model allows to implement the SUSY-based mechanism 
of [14, 15] that ensures emergence of LI at low energies within the SM sector. SUSY turned 
out to be more powerful in this respect than one could a priori expect: it leads to emergent 
low-energy LI even within the aether sector itself. We have found that the dynamics of the 
super-aether is richer than of its non-SUSY counterpart. In particular, the model possesses a 
family of inequivalent vacua exhibiting different symmetry breaking patterns. Remarkably, 
the model tolerates a certain breaking of spatial isotropy while remaining stable and ghost 
free. Finally, we have analyzed the effects of SUSY breaking. 

In this paper we have restricted to the case of global SUSY. An important next step is 
the coupling of the super-aether model to supergravity. This task appears straightforward 
but presents a non-trivial technical exercise. Leaving it for the future, let us mention the 
following subtlety. It is impossible to embed the model of the present paper in the framework 
of the minimal SUGRA [25, 26]. The reason is that in this framework, due to non- vanishing 
super curvature, the covariant spinor derivatives do not anticommute when acting on a vector 
superfield, 

{Va.V^^O. (47) 

This prevents the definition of the chiral aether vector [27]. To overcome this obstruction 
one has to resort to the non-minimal formulation of SUGRA [28, 29] where it is possible to 
construct improved spinor derivatives, that do anticommute [30]. 

It would be interesting to generalize the construction of this paper to other theories with 
dynamical preferred frame. To illustrate possible problems let us consider the khrono- metric 
model of [9]. This model arises as the low-energy limit of the healthy Hofava gravity [6, 7] 
and is closely related to the Einstein-aether theory [31]. Violation of LI in this model is 

11 Strictly speaking, the SUSY breaking will also affect the relations (33) through higher-dimensional 
operators. However, these corrections are suppressed by powers of the ratio of the SUSY breaking scale to 
Ma,. They are small provided there is a hierarchy between these scales. 
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described by a real scalar field ip - the khronon - that has non- vanishing time-like gradient. 
One imposes invariance under reparameterizations of <p, 



p ^ <p((p) , (48) 

which implies that <p enters the Lagrangian only through the unit time-like vector 

d m (p 
v / -d n (pd n (p 

In terms of this vector the action has the same form (3) as the aether action 12 . Note, however, 
that in terms of the fundamental field ip the action contains higher derivatives leading to 
fourth-order equations of motion. As discussed in [9], this does not lead to inconsistencies. 
Due to the special structure of the fourth-order operator it is always possible to choose a 
coordinate frame in such a way that the equations contain only two time-derivatives, the extra 
two derivatives being purely space-like. The preferred frame is determined by the khronon 
field itself. Namely, the time coordinate must be chosen to coincide with the background 
value of the field, 

t — ^background • (50) 

We stress that the possibility to make this choice is crucial for the consistency of the model. 
As we now discuss, this property is lost when we attempt to construct a SUSY extension of 
the model, and higher time derivatives persist in the equations of motion. 

By the analogy with the aether case, a natural choice of the khronon superfield is a chiral 
scalar, 

$ = <p(y) + V29 a Mv) + 2 F(y) . (51) 

Importantly, this requires that the khronon field ip be complexified. From this scalar one 
constructs a vector superfield 

U m = -=M=. (52) 

The lowest components of U m and $ are related by Eq. (49). Note that U m is automatically 
chiral and satisfies the unit-norm constraint (22). The Lagrangian of the theory has the form 
(24) with the only difference that now there is no need for a Lagrange multiplier. 
The khronon vacuum corresponding to (35) has the form 

fvac — t cos a + iz sin a . (53) 



12 For the hypersurface-orthogonal vector (49) the four terms in the aether action are not independent. 
One of them can be eliminated in favor of the rest. 
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Note that this vacuum preserves SUSY despite of the fact that the SUSY variation of the 
fermionic component is non-zero, 

Sip = iV2a m ( d m (p vac = V2(ia° cos a - a 3 sin a)( , (54) 

where ( is the parameter of the SUSY transformation. The reason is the large symmetry 
group of the theory. The expression (52) is invariant not only under arbitrary reparameter- 
izations of $, but also under all transformations of the form 

0°) (55) 

that involve arbitrary dependence on the holomorphic Grassmann coordinates. An appro- 
priate transformation of this form compensates the variation (54). 

Consider now small perturbations on top of the vacuum (53). For our argument it is 
sufficient to concentrate on the perturbations of the khronon, leaving aside the fermionic 
and auxiliary fields. We write 

V = Vvac + X , (56) 

and expand the action to quadratic order in x- We can use the result (38) for the aether: 
it remains to express the aether perturbations V\^i v \\ m terms of \- A straightforward 
calculation yields 

v a = d a x , a = 1,2, (57a) 
v\\ = cos a d 3 x — ? sin a x • (57b) 

If a = we recover the quadratic Lagrangian of [9] for the (complexified) khronon pertur- 
bations containing only first time derivatives of x- However, in the general case a ^ the 
action involves second time derivatives x an d results in the fourth-order equations of motion. 
This leads to very fast instabilities that are clearly unacceptable. Note that the problem 
cannot be eliminated even by breaking SUSY or simply by choosing the background with 
a = 0, since the higher time derivatives will inevitably appear in interaction terms. The ori- 
gin of this behavior can be traced back to the fact that the khronon is complex and therefore 
there is no real coordinate system in which it coincides with the would-be preferred time. 
Thus we conclude that new ideas are required to supersymmetrize the khronon model and, 
eventually, its UV completion - Hofava gravity. 

Another issue related to the problem of constructing a SUSY Hofava gravity, but also 
interesting on its own right, is the extension of the notion of SUSY to theories invariant 
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under anisotropic scaling of space and time coordinates. Recently, it has been demonstrated 
[24] that the requirements of the anisotropic scaling and gauge invariance are incompatible 
with the non-relativistic SUSY algebra (12). Namely, it is impossible to construct any 
gauge theory 13 with anisotropic scaling invariant under (12). As a way out one could try 
to deform the SUSY algebra by allowing non-linear functions of the spatial momentum P a 
appear on the r.h.s. of the commutation relations [23]. This would make SUSY compatible 
with the anisotropic scaling. However, as pointed out in [24], construction of interacting 
models invariant under the deformed algebra may present a non-trivial task because the 
SUSY transformations now contain higher spatial derivatives. In particular, one cannot 
use the superfield formalism to build invariant Lagrangians. Indeed, any attempt to realize 
the supercharges on the superspace will involve higher-order differential operators. These 
supercharges will not satisfy the Leibniz rule, implying that the product of two superfields 
does not transform as a superfield and thus the superfield formalism becomes useless. At 
present only free theories invariant under the deformed SUSY have been constructed [24]. 
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